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Abstract 

Integrable systems are usually given in terms of functions of con- 
tinuous variables (on R), in terms of functions of discrete variables 
(on Z), and recently in terms of functions of q- variables (on We 
formulate the Gel'fand-Dikii (GD) formalism on time scales by us- 
ing the delta differentiation operator and find more general integrable 
nonlinear evolutionary equations. In particular they yield integrable 
equations over integers (difference equations) and over g-numbers (q- 
difference equations) . We formulate the GD formalism also in terms of 
shift operators for all regular-discrete time scales. We give a method 
allowing to construct the recursion operators for integrable systems 
on time scales. Finally, we give a trace formula on time scales and 
then construct infinitely many conserved quantities (Casimirs) of the 
integrable systems on time scales. 







1. Introduction 



Integrable systems are well studied and well understood in 1 + 1 dimensions 
[l]-[3]. Here one of the dimensions denotes the time (evolution) variable and 
the other one denotes the space variable which is usually taken as continuous. 
There are also important examples where this variable takes values in Z, i.e., 
integer values. In both cases the Gel'fand- Dikii (GD) approach is quite 
effective. One can generate hierarchies of integrable evolution equations, 
both on K. and on Z (see [3] for GD applications and related references). 
In addition one can construct the conserved quantities, Hamilton operators, 
and recursion operators. Investigation of integrable systems on g-discrete 
intervals started in [4]- [6] . They considered GD formalism on K q and found 
g-integrable hierarchies including the g-KdV equation. 

In this work we extend the Gel'fand- Dikii approach to time scales where 
R, Z, and IK^ are special cases. In the next section we give a brief review of 
time scales calculus. See the references [7]- [13] for a more detailed review of 
the subject . In GD formalism, in obtaining integrable systems the essential 
tools are the differential and shift operators and their inverses. For extending 
the GD formulation to time scales we give necessary means to construct in the 
sequel the algebra of pseudo A— differential operators and the algebra of shift 
operators. In Section 3 we assume A-differential Lax operators and derive 
the A-Burgers hierarchy with its recursion operator. We present special cases 
of the Burgers equation for T = hZ and T = K q . In Section 4, we consider 
the regular time scales where the inverse of jump operators can be defined. 
Here we assume a pseudo A-differential algebra and give the corresponding 
GD formulation. As an example we present a A-KdV hierarchy. We first 
find n = 1 member of the hierarchy and write out it explicitly for T = 
R, Z, K q and for T = (— oo, 0) U K q . Then we give the n = 3 member 
and call it as the A-KdV system. We call it A-KdV equation, because 
the corresponding Lax operator is a second order A-differential operator. 
It involves two fields u and v, but the second field v can be expressed in 
terms of the first filed u. When T = R, this system reduces to the standard 
KdV equation. In Section 5, we consider the regular- discrete time scales and 
introduce the algebra of shift operators on them and give the corresponding 
GD formulation for all such time scales. Here several examples are presented. 
We first generalize the examples of discrete systems on Z given in [3] (one 
field, two fields, and four fields examples in [3]) to arbitrary discrete time 
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scales. In all these examples when T = Z we get the discrete evolutions 
given in [3] . We construct the recursion operators of these systems on time 
scales. We generalize the Frenkel's KdV system [4] introduced on K q to 
arbitrary discrete time scales and we construct its recursion operator. In 
this section, we finally give an example of the KP hierarchy on discrete 
time scales. In Section 6, we extend the standard way of constructing the 
conserved quantities of integrable systems to time scales by introducing a 
trace form on the algebra of A-pseudo differential operators. The trace form 
introduced in this section reduces, in particular cases, to the standard trace 
forms on R and Z. In Appendix we give the recursion operators of two four 
fields systems introduced in Section 5. We end up with a conclusion. 

2. Time Scale Calculus 

The time scale calculus is developed mainly to unify differential, difference, 
and q— calculus. A time scale (T) is an arbitrary nonempty closed subset of 
the real numbers. The calculus of time scales was initiated by Aulbach and 
Hilger [7], [8] in order to create a theory that can unify and extend discrete 
and continuous analysis. The real numbers (M), the integers (Z), the natural 
numbers (N), the non-negative integers (No), the /i-numbers (/iZ = {hk : 
k G Z}, where h > is a fixed real number), and the q— numbers {K q = 
q z U {0} = {q k : k G Z} U {0}, where q > 1 is a fixed real number) are 
examples of time scales, as are [0, 1] U [2,3], [0, 1] U N, and the Cantor set, 
where [0,1] and [2,3] are real number intervals. In [7], [8] Aulbach and 
Hilger introduced also dynamic equations on time scales in order to unify 
and extend the theory of ordinary differential equations, difference equations, 
and quantum equations [9] (/i-difference and q -difference equations based on 
/i-calculus and g-calculus,respectively). For a general introduction to the 
calculus on time scales we refer the reader to the textbooks by Bohner and 
Peterson [10], [11]. Here we give only those notions and facts connected to 
time scales which we need for our purpose in this paper. 

Any time scale T is a complete metric space with the metric (distance) 
d(x,y) = \x — y\ for x,y G T. Consequently, according to the well-known 
theory of general metric spaces, we have for T the fundamental concepts 
such as open balls (intervals), neighborhood of points, open sets, closed sets, 
compact sets, and so on. In particular, for a given number r > 0, the r- 
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neighborhood U r (x) of a given point x G T is the set of all points y G T 
such that d(x, y) < r. By a neighborhood of a point x G T is meant an 
arbitrary set in T containing an r-neighborhood of the point x. Also we 
have for functions / : T — > R the concepts of the limit, continuity, and 
properties of continuous functions on general complete metric spaces (note 
that, in particular, any function / : Z — > R is continuous at each point of Z). 
The main task is to introduce and investigate the concept of derivative for 
functions / : T — > R. This proves to be possible due to the special structure 
of the metric space T. In the definition of derivative, the so-called forward 
and backward jump operators play special and important roles. 

Definition 1. For x G T we define the forward jump operator a : T — > T by 

a(x) =inf {y GT : y > x}, (1) 

while the backward jump operator p : T — > T is defined by 

p{x) = sup {y G T : y < x}. (2) 



In this definition we put in addition <r(maxT) = maxT if there exists a finite 
maxT, and p(minT) = minT if there exists a finite minT. Obviously both 
a(x) and p(x) are in T when x G T. This is because of our assumption that 
T is a closed subset of R. 

Let x G T. If o~{x) > x, we say that x is right- scattered, while if p(x) < x 
we say that x is left- scattered. Also, if x < maxT and o~(x) = x, then x is 
called right-dense, and if x > minT and p(x) = x, then x is called left-dense. 
Points that are right-scattered and left-scattered at the same time are called 
isolated. Finally, the graininess functions p, v : T — > [0, oo) are defined by 

p(x) = a(x) — x, and v(x) = x — p(x) for all iGl (3) 

Example 1. If T = R, then a(x) = p(x) = x and p(x) = v(x) = 0. If 
T = hZ, then a(x) = x + h, p(x) = x — h, and p(x) = v{x) = h. On the 
other hand, if T = K q then we have 

a(x)—qx, p(x) = q~ l x, p(x) — (q — l)x, and v{x) = (1 — q' 1 ) x. (4) 
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Let T K denotes Hilger's above truncated set consisting of T except for 
a possible left-scattered maximal point. Similarly, T K denotes the below 
truncated set obtained from T by deleting a possible right-scattered minimal 
point. 

Definition 2. Let f : T — > K. be a function and x G T K . Then the delta 
derivative of f at the point x is defined to be the number f A (x) (provided it 
exists) with the property that for each e > there exists a neighborhood U of 
x in T such that 

\f(a(x)) - f(y) - f A (x)[a(x) - y]\ < s\a(x) - y\, (5) 

for all y G U . 

Remark 1. If x G T \ T K , then f A (x) is not uniquely defined, since for 
such a point x, small neighborhoods U of x consist only of x and besides we 
have a(x) = x. Therefore (5) holds for an arbitrary number f A (x). This is 
a reason why we omit a maximal left-scattered point. 

We have the following: (i) If / is delta different iable at x, then / is con- 
tinuous at x. (ii) If / is continuous at x and x is right-scattered, then / is 
delta differentiable at x with 

(iii) If x is right-dense, then / is delta differentiable at x iff the limit 

Hm /W-/M) 

y->x x — y 

exists as a finite number. In this case f A (x) is equal to this limit, (iv) If / 
is delta differentiable at x, then 

f(a(x)) = f(x) +f ,(x)f A (x). (8) 



Definition 3. If x G T K , then we define the nabla derivative of f : T — » 1R at 

x to be the number f v (x) (provided it exists) with the property that for each 
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e > there is a neighborhood U of x inT such that 

\f(p(x)) - f(y) - f{x)[p{x) - y]\ < e\p{x) - y\, (9) 

for all y G U . 

We have the following: (i) If / is nabla differentiable at x, then / is contin- 
uous at x. (ii) If / is continuous at x and x is left-scattered, then / is nabla 
differentiable at x with 

A*) = f{x) ~£f {x)) - (io) 

(hi) If x is left-dense, then / is nabla differentiable at x if and only if the 
limit 

hm f{x) - fiy) (11) 
y^x x — y 

exists as a finite number. In this case f v (x) is equal to this limit, (iv) If / 
is nabla differentiable at x, then 

f(p(x))=f(x)-v(x)f v (x). (12) 

Example 2. If T = R, then f A (x) = f v (x) = f'(x), the ordinary derivative 
of / at x. If T = hZ, then 

fA{x) ^ f(x + h)^-f(x) ^ r{x) = f(*)-f(*-h) m (13) 
If T = K q , then 

^-Sfc™ wd n ^/M^l£), (14 , 



for all x 7^ 0, and 



/ A (0) = / v (0) = lim /(0) (15) 

y->0 y 



provided that this limit exists. 
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Among the important properties of of the delta differentiation on T we have 
the Leibnitz rule: If /, g : T — > R are delta differential) le functions at x G T K , 
then so is their product fg and 

(fg) A (x) = f A (x)g(x) + f(a(x)g A (x), (16) 
= f(x)g A (x) + f A (x)g(a(x)). (17) 

Also, if /, g : T — > R are nabla differentiable functions at rr G T K , then so is 
their product fg and 

(fg) v (x) = f V (x)g(x) + f(p(x)g v (x), (18) 
= f(x)g^(x) + r(x)g(p(x)). (19) 

In the next proposition we give a relationship between the delta and nabla 
derivatives (see [12]). 

Proposition 4. (i) Assume that f : T — > R zs c?e/to differentiable on T K . 
T/ien / nafr/a differentiable at x and 

f(x) = f A (p(x)), (20) 

for x G T K swc/i £/ia£ o"(p(x)) = x. 7/j m addition, f A is continuous on T K , 
then f is nabla differentiable at x and (20) holds for any xeT s . 
(ii) Assume that f : T — > R na&Za differentiable on T K . Then f is delta 
differentiable at x and 

/ A (x) =/>(*)), (21) 

for x G T K snc/i i/iat p(a(a;)) = x. 7/j m addition, / v zs continuous on T K , 
then f is delta differentiable at x and (21) holds for any x G T K . 

Now we introduce the concept of integral for functions / : T — > R. 

Definition 5. ^4 function F : T — > R called a A-antiderivative of f : T — > 
R provided F A (x) = f(x) holds for all x inT K . Then we define the A-integral 
from a to b of f by 

j f(x)Ax = F(b)-F(a) for all a, b G T. (22) 

J a 



Definition 6. A function $ : T — > R called a V ' -antiderivative of f : T — > 
R provided $ v (x) = /(x) /ioWs /or a// x m T K . T/ien we define the V '-integral 
from a to b of f by 

j f( x )Vx = $(b)-$(a) for all a,beT. (23) 

If a, fe G T with a < b we define the closed interval [a, 6] in T by 

[a,b] = {x G T : a < x < b}. (24) 

Open and half-open intervals etc. are defined accordingly. Below all our 
intervals will be time scale intervals 

Example 3. Let a, b G T with a < b. Then we have the following: 

(i) If / : T -> R then 

/ f(x)Ax= [ f{x)Vx= f f(x)dx, (25) 

J a J a J a 

where the integral on the right-hand side is the ordinary integral. 

(ii) If [a, b] consists of only isolated points, then 

/ f(x)Ax = vi x )fi x ) and / f(x)Vx= »i x )fi x )- ( 26 ) 

Ja xe[a,b) Ja xe{a,b] 

In particular, if T = Z, then 

f f(x)Ax = Y^f(k) and f/(x)Vx= ]T f(k). (27) 

Ja k=a Ja k=a+l 

If T = hZ, then 

/ f(x)Ax = h f( x ) and / f(x)Vx = h f( x ) ( 28 ) 

^ a xe[a,b) ^ a xe(a,b] 
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and if T = K ? , then 

/b pb 
f(x)Ax = (l-q) ^2 x f( x ) and / /(i)Vi=(l-g' 1 ) x f( x )- 

xe[a,b) ^ a xe(a,b] 

(29) 

The following relationship between the delta and nabla integrals follows 
from Definitions 5 and 6 by using Proposition 4. 

Proposition 7. // the function f : T — > R zs continuous, then for all a, b G T 

a < b we have 

f f(x)Ax= [ f(p(x))Vx and [ f{x)Vx= f f(a(x))Ax. (30) 



Indeed, if F : T — > R is a A-antiderivative for /, then F A (:r) = /(x) for all 
x G T K , and by Proposition 4 we have f(p(x)) = F A (p(x)) = F v (x) for all 
x G T K , so that F is a V-antiderivative for f(p(x)). Therefore 

f f(p(x))Vx = F(b) - F(a) = f f(x)Ax. (31) 

J a J a 

From (16)-(21) and (30) we have the following integration by parts formu- 
las: If the functions /, g : T — > R are delta and nabla differentiable with 
continuous derivatives, then 



I" f A (x)g(x)Ax = 

J a 


f(x)9(x)\ b a 


- f f(a(x))g A (x)Ax, 

J a 


(32) 


f f{x)g{x)Vx = 

J a 


f( x )g( x )\ b a 


~ f ' f(p(x))gV( x )Vx, 

J a 


(33) 


f f A (x)g(x)Ax = 


f( x )g( x )\ b a 


- f f(x)g v (x)Vx, 

J a 


(34) 


f f{x)g{x)Vx = 

J a 


f( x )g( x )\ b a 


- f f(x)g A (x)Ax. 

J a 


(35) 


For more general treatment of the delta 


integral on time scales (Riemann 



and Lebesgue delta integrals on time scales) see [13] and Chapter 5 of [11]. 
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3. Burgers Equation on Time Scales 



The Gel'fand-Dikii approach is very effective in studying the symmetries, 
bi-Hamiltonian formulation, and in constructing the recursion operators of 
integrable nonlinear partial differential equations. In this approach one takes 
the Lax operator L in an algebra like a differential or pseudo differential 
algebra, a matrix algebra, a polynomial algebra, or the Moyal algebra. In 
this section we take L in the algebra of delta differential operators. 

Let T be a time scale. We say that a function / : T — > M. is A-smooth if 
it is infinitely A-differentiable (and hence infinitely V-differentiable). By 
A we denote the delta differentiation operator which assigns to each A- 
differentiable function / : T — > R its delta derivative A(/) defined by 

[A{f)]{x) = f A {x) for xeT K . (36) 

The shift operator E is defined by the formula 

(Ef)(x)=f(a(x)) (37) 

for x G T, where a : T — > T is the forward jump operator. It is convenient, 
in the operator relations to denote the delta differentiation operator by 5 
rather than by A. For example, 5f will denote the composition (product) of 
the delta differentiation operator 8 and the operator of multiplication by the 
function /. According to formula (16)) we have 

Sf = f A + E(f)5. (38) 
Consider the N-ih order ^-differential operator given by 

L = a N 5 N + oat-i 5 JV ~ 1 + • • • + ai 5 + a Q , (39) 

where the coefficients (i — 0, 1, • • • N) are some A-smooth functions of the 
variable i6l These functions are assumed to depend also on a continuous 
variable i G 1, however, we will not (for simplicity) indicate explicitly the 
dependence on t. 

Proposition 8. Let L be given as in (39) and A n = (L n ) >0 be the operator 
L n missing the 5° term. Then the Lax equation 
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^ = [A n ,L]=A n L-LA n (40) 
dt n 

for n = 1, 2, • • • produces a consistent hierarchy of coupled nonlinear evolu- 
tionary equations. 

Example 4. Burgers equation on time scale: Let L — v 8 + u, where u 
and v are functions of x G T and t G R. Then for an appropriate operator A 
the Lax equation 

§ = [A,L] (41) 

defines a system of two differential equations for the functions u and v. We 
find the operator A by using the Gelfand-Dikii formalism. Let us start with 
the second power of L and assume A = (L 2 ) >0 , where 

L 2 = v E(v) 5 2 + [vv A + v E(u) + uv] 5 + v u A + u 2 . (42) 

We can assume A = —(L 2 ) (the part of — L 2 without the 5 terms). With 
this choice, (41) gives 

| = ^(vu A + u 2 ) A , (43) 

^ = v (vu A + u 2 ) A , (44) 
at 

where fi(x) = o~(x) — x for x G T. 

Equations (43) and (44) given above are not independent of each other. It is 
easy to see that v = /iu + A, where A is an arbitrary real function depending 
only on x G T. Then these two equations reduce to a single equation, a 
Burgers equation on time scales, 

— = ( fiU + \)[u 2 + (piu + \)u A ] A . (45) 
at 

Let us present some special cases: (i) When T = R then /j, — and 5 = D, 
the usual differentiation. Hence we can let A = 1 and (45) reduces to the 
standard Burgers equation on R.(ii) When T = KL then //(to) = h and 
f A (m) = \[f(m + h) — /(to)] for any /. Then taking A = in (45) we find 
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du(m) 

— - — = u(m) u(m + h) [u(m + 2h) — u(m)], (46) 
at 

where m G KL. The evolution equation given above in (46) represents a 
difference version of the Burgers equation, (iii) Let T = q z , where g ^ 1 
and q > 0. Then we have p(x) — (q — \)x and f A (x) = an d taking 

A = we get from (45) 

^ u ( x ^ — u(x)u(qx)[u(q 2 x) — u(x)]. (47) 

Taking A n = — {L n )o with L given as in Example 4 we get a hierarchy of 
evolution equations (Burgers hierarchy on time scales) from 

£ = -[(L») , ,,* + „] (48) 

for all n — 1, 2, 3, • • •. Since (L n ) is a scalar function, letting (L n ) = p n we 
obtain 



^ = W( Pn ) A , (49) 
^ = -(Pn) A , (50) 

where the first three p n are given by 

Pi = u, (51) 
p 2 = VM A + M 2 ? (52) 

p 3 = w £(t;) M AA + [vv A + vE(u) + ™]m a + {vu A + u 2 )u. (53) 
The above hierarchy reduces to a single evolution equation with v — pu + A: 

^- = (pu + X)(p n ) A , n=l,2,... (54) 

where p„ is equal to p n with t> = pu + A. When T is a regular-discrete time 
scale, the first three p n are given for A = by 
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Pi = u, (55) 
p 2 = uE(u), (56) 
p 3 = uE(u)E 2 (u). (57) 



It is possible to construct the recursion operator 1Z by using the Lax 
representation [14]- [16]. The hierarchy satisfies a recursion relation like 

dL dL r . 

L — + [R n ,L], n = l,2,.-. (58) 



where i? n is the remainder operator which has the same degree as the Lax 
operator L. We shall construct this operator for the Burgers equation with 
A = on regular-discrete time scales. Choosing R n = a n 5 we get (by 
choosing v(x) = ji{x) u(x)) 

n = uE+[u {E{u) - u)\ (E - I)" 1 (59) 

E[u) 

One can generate the hierarchy (54) by application of the recursion operator 
1Z to the lowest order symmetry U\ — u (E(u) — u): 

^ = K n - x u u n = l,2,.... (60) 



4. Algebra of Pseudo Delta Differential Oper- 
ators on Regular Time Scales 

Let us define the notion of regular time scales. 

Definition 9. We say that a time scale T is regular if the following two 
conditions are satisfied simultaneously : 

(i) a(p(x)) = x for all x e T and (61) 

(ii) p(a(x)) = x for all x e T, (62) 
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where a and p denote the forward and backward jump operators, respectively. 

From (61) it follows that the operator a : T — > T is "onto" while (62) im- 
plies that a is "one-to-one". Therefore a is invertible and o~ x = p. Similarly, 
the operator p : T — > T is invertible and = a if T is regular. 

Let us set = minT if there exists a finite mini, and set = — oo 
otherwise. Also set x* = maxT if there exists a finite maxT, and x* = oo 
otherwise. It is not difficult to see that the following statement holds. 

Proposition 10. A time scale T is regular if and only if the following two 
conditions hold: 

(i) The point x* = T is right-dense and the point x* = maxT is left-dense. 

(ii) Each point o/T\{a;*, x*} is either two-sided dense or two-sided scattered. 

In particular, R, /iZ, and K g are regular time scales, as are [0, 1], [—1, 0]U{^ : 
k E N} U : k e N} U [1, 2], and (-oo, 0] U {± : k E N} U {2k : k E N}, 
where [—1, 0], [0, 1], [1,2], (— oo, 0] are real line intervals. 

If / : T — > R is a function we define the functions f a : T — > R and 

: T -> R by 

/ CT ( X ) = /((j(x)) and / p (x) = for all x E T. (63) 

Defining the shift operator E by the formula E f = f a we have 

(£/)(a;) = f a (x) = f(a(x)) for all x E T. (64) 

The inverse E~ x exists only in case of regular time scales and is defined by 

(E- 1 f){x) = f(a-\x)) = f(p(x)) for all x E T. (65) 

In the operator relations, for convenience, we will denote the shift operator by 
£ rather than by E. For example, £ f will denote the composition (product) 
of the shift operator £ and the operator of multiplication by the function /. 
Obviously, for any integer m E Z, we have 

£ m f = (E m f)£ m . (66) 

Remember that 5 denotes the delta differentiation operator acting in the 
operator relations by 5f = f A + E(f)S. The following proposition is an 
immediate consequence of the the formulas (8), (16). 



13 



Proposition 11. The operator formulas 

£ — I + /i5, and (67) 
5f = f A + E(f)5 (68) 

hold, where the function fi : T — > R defined by /i(x) = a(x) — x for all 
x£T, and / denotes the identity operator. 

In this section we will assume that all our considered functions from T to 
M are A-smooth and tend to zero sufficiently rapidly together with their A- 
derivatives as x goes to x* or x*, where = minT if there exists a finite 
minT and x* = — oo otherwise, x* = maxT if there exists a finite maxT 
and x* = oo otherwise. The inverse operator exists on such functions. If 
g : T — > T is such a function, then 

[A" 1 (^)]W= /%(y)Ay. (69) 

Proposition 12. Lei / : T -> R 6e a A-smooth function such that f and all 
its A-derivatives vanish rapidly at and x* . Then the operator 5~ l f being 
the composition (product) of 5~ l and f has the form of the formal series in 
powers of 5~ l 

(T 1 / = a <r 1 + ai 5~ 2 + ■■■ (70) 
where a = E~ l f, and a k = (-l) k (E' 1 f) yk for k = 1, 2, • • •. 
Proof. Multiplying (68) on the left and right by <5 -1 we obtain 

S~ 1 E(f) = fS' 1 -r 1 f A S-\ (71) 
Replacing here / by E^ 1 f we get 

r 1 / = (E- 1 f) r 1 - 5-\E- 1 f) A 8- 1 . (72) 

Further, applying this rule to the function (E~ 1 f) A and taking into account 
that by Proposition 4(i) 

E- 1 (E- 1 f) A = (E-y , (73) 

we find 
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Substituting this into the second term on the right-hand side of (72) we 
obtain 

r 1 / = z)*- 1 - (£~7) v r 2 + r 1 ((£"7) A ) V r 2 . (75) 

Continuing this procedure repeatedly we arrive at the statement of the propo- 
sition. 

Definition 13. By A we denote the algebra of pseudo delta differential 
operators. Any operator K e A of order k has the form 

k 

k = a ^ e ( 76 ) 



where ag's are A-smooth functions of x e T. For K given by (76) we will 
use the following notations: 

k -l 
K> = a e 5 e , and K <0 = ^ a e S e . (77) 
e=o -oo 



As an example we let 

L — a N 5 N + a N -\ 5 N ~ 1 H ^a 1 5 + a (78) 

where ai(i = 0, 1, • • • , N) are some A-smooth functions on T. Then we have 

Proposition 16. Let L be given in (78). For each fixed N the Lax equation 

^ = [A n ,L], A n = (L%) >0 , (79) 
at n 

for n = 1, 2, • • • not divisible by N, produces a (consistent) hierarchy of evo- 
lution equations (a KdV hierarchy on time scales) . 

Proof. Since (L^)> = L% — (L^) <0 , we get 
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^ = [(L*)> ,L] = -[L*) <0 ,L]. (80) 

Evidently the commutator [(L n//Ar )> , L] involves only nonnegative powers of 
5, while the commutator [(L n ^ N ) < Q, L] has the form Ylf=~oo bj & ■ Therefore, 
we get by (80) that, for all n not divisible by N , (79) produces nontrivial 
consistent N + 1-number of evolutionary coupled A-differential equations 
for aj, i = 0,1, • • • , N. Note that turns out to be a fixed (i.e., time 
independent) function of x. 

Example 5. A KdV hierarchy on time scales. Let 

L = 5 2 + v5 + u, (81) 
where u and v are A-smooth functions. It is straightforward to find that 

L 1 / 2 = 5 + a + a 1 5~ 1 + a 2 5- 2 + ■■■ (82) 

where 

E(a ) + a = v, (83) 
J B(a 1 ) + a 1 + (a ) A + (a ) 2 = M, (84) 
E(a 2 ) + a 2 + a l E- 1 (a ) + {a l )^ = 0. (85) 

Choosing n — 1, 3, • • • we get the members of the KdV hierarchy. 
(1). Let n — 1. Then Lax equation (79) becomes 

and gives coupled equations for u and v 

f = u A -v(a ) A -(a ) AA , (87) 

— = v A + E(u) — u — v [E(a ) — a ] 

-E{a A ) - E(a ) A , 
= ^(u A -v(a ) A -(a ) AA ). (88) 



16 



Comparing the above equations we get 

dv du , . 

*-"w = - (89) 

and therefore 

v = fiu + X, (90) 

where A is an arbitrary real function depending only on x G T. Thus, two 
equations (87) and (88) reduce to the following single equation 

*i = M A _ {fxu + A) (afl) A _ (ajj) AA j (Q1) 

where a is expressed, according to (83), from 

E(a ) + a = fiu + A. (92) 
If we take A = 0, then (91) and (92) become 

^ = M A -^(a ) A -M AA , (93) 
E(a ) +a = fiu. (94) 

We shall now give a , for illustration, for particular cases of T: 

(i) In the case T = lwe have fx — and (94) gives a = and (93) becomes 

du du _ . 

* = (95) 

which is a linear equation explicitly solvable: 

u(x,t) = + (96) 

where ip is an arbitrary differentiable function. 

(ii) In the case T = Z we have fj, — 1 and (94) is satisfied by 

a (n) = - ^ (-l)" +fe M(A;), nGZ (97) 



k=— oo 
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and therefore the equation (93) becomes 

d ^ = - u \n) + 2u{n) + 2{-l) n [2 + u{n)} J2(-l) k u (k), (98) 

k=— oo 

for n G Z. 

(iii) In the case T = K q we have /i(x) = (q — l)x and (94) is satisfied by 
ao(0) = and 

a (x) = -(q-l) (-iy° 9qixy) yu{y) (99) 

for x G Kg and x ^ 0. Substituting (99) into (93) we can get an evolution 
equation for u. 

(iv) Let T = (-oo, 0) U K q = (-oo, 0] U q z . In this case fi{x) = if 
x G (— oo, 0] and /i(x) — (q — l)x if x G g z . The equation (94) is satisfied by 
the function a given by 



^ " I -(? - 1)2^-1,] (-1)^) yufc,) x G g z 

Therefore (93) will yield an evolution equation coinciding on (— oo, 0] and q z 
with the evolution equations described in the examples (i) and (iii), respec- 
tively. Now an essential complementary point is that the solution u must 
satisfy at x = the smoothness conditions 

u(0") = w(0 + ), u'(O-) = u A (0 + ). (101) 
(2). Letting n — 3, first we get 

^3/2 = ^3_|_^^2_|_^^_|_ r _|_ ^ erms w j^] 1 negative powers of S) (102) 



where 



p = a + E(v), (103) 
q = V A + S(u) + a u + on, (104) 
r = m a + a>o u + ai E~ x (v ) + a 2 , (105) 
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and the Lax equation 



* f + | = [( tf/> )i0 , £] , (106 ) 



gives the coupled equations for u and v 

d u AAA , AA , A AA A / in7 \ 

— = u + pu + qu — r —vr , (107) 

(JjL 

^ = v AAA + E(u AA ) + (E(u A )) A + E(u) AA +p[v AA + E(u A ) 
at 

+E(u) A ] + q(v A + E{u) -u)+rv- q AA - E(r A ) 

-E(r) A -vq A - vE(r). (108) 



As in the first member of the hierarchy (n — 1 case), the above A-KdV 
equations reduce to a single equation for the function u. Below in Corollary 
23 we found that v = /i(x)u + X(x). Letting A = constant we get 

d u AAA , AA , A AA A /mn\ 

— — u+pu+qu—r—vr. (109) 

(JjL 

It is possible to write the above equation more explicitly in terms of u for 
T = R, T = Z, and for T = K q but they are quite lengthy. For the discrete 
case we give a KdV hierarchy in Example 8, next section. 



5. Shift Lax Operators on Regular-Discrete 
Time Scales 

Let T be a time scale. Let us set x* = minT if there exists a finite minT 
and x* = —oo otherwise. Also set x* = maxT if there exists a finite maxT 
and x* = oo otherwise. We will briefly write x* = minT and x* = maxT. 

Definition 17. We say that a time scale T is regular- discrete if the following 
two conditions are satisfied: 

(i) The point Ob ^ %S right-dense and the point x* is left-dense. 

(ii) Each point of T\{x*, x*} is two-sided scattered (isolated). 

The shift operator E is defined on functions / : T — > R by the formula 
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(Ef)(x) = f(a(x)), for 16T (110) 

where a : T — > T is the forward jump operator. 

In this section we deal only with regular-discrete time scales T. For such 
time scales T we have 

li(x) — <j(x) — x ^ 0, for all x £ T\{x*,x*} (HI) 

and, therefore, on functions given on T\{x*,x*} we have the operator rela- 
tionship 

6=-(S-l). (112) 

All our functions will be assumed to be defined on T\{x*,x*} and tends to 
zero sufficiently rapidly as x goes to x* or x*. 

This shift operator £, should be quite useful in the application of the 
Gel'fand-Dikii formalism. The reason is that for any integer m we have the 
simple product rule 

£ m u = (E m u) £ m . (113) 

Hence, for regular-discrete time scales, we can define an algebra of £ opera- 
tors. 

Definition 18. An algebra , A e7 of £ operators satisfying the operator equa- 
tion (113) is defined as follows: Any operator K in A e with degree k is of the 
form 

k 

K = ^az£ l (114) 

— oo 

where ae are functions of x e T that depend also ont<E~R. 

Hence we can form Lax operators in A e , and produce integrable equations 
on regular-discrete time scales. Following [3] we obtain two classes of Lax 
representations. 

Proposition 19. The Lax equation 

% = [(L e ) >k ,L], A; = 0,1 (115) 
dtt 
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produces consistent hierarchy of equations for £ = 1, 2, • • • with the following 
suitable Lax operators 



L = S a+n + u a+n ^ 1 S a+n - 1 + --- + u a S a , (116) 
L = Va+n S a+n + v a+n ^S a+n - 1 + --- + v a+l S a+l + S a , (117) 

for k = and k = 1, respectively. Here u,i and Vi are functions defined on T 
and the integer a is restricted to satisfy the inequality — n < a < — 1. 

Remark. Lax operators above and the following examples are given on any 
regular-discrete time scale T (we can take in particular T = Z or K 9 ). This 
means that for any function u on such a time scale E{u) = u(a(x)) where 
a is the jump operator defined in the second section. Hence our examples 
and results should be considered as more general than those considered in 
[3]. In the case of Ref.[3] time scale is just the integers (T = Z) where 
E(u(n)) = u(n + l) 

Example 6. Two field equations. Let k — 0, a — — 1 and 

L = u_ 1 S' 1 + u + S = v£- 1 + u + S. (118) 



Then we find 



£ = 1 ^ = v{u-E-\u)), (119) 
IT, = E{V) - V > (120) 



dv 

£ = 2 — = ^v + E^v-vE^^-vE' 1 ^ 2 ), (121) 



dt 2 
du 

df 2 



uE(v) + E(u)E(v) -vE-\u) -uv, (122) 



dv 

£ = 3 — = uv 2 + u*v-vE~ l (v)E- 2 (u)-2vE- l (u)E- l (v) 



dt 



3 
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vE-\u 3 ) + 2uvE(v) - v 2 E- 1 {u) + vE(u)E(v), (123) 

^ = E(v)[u 2 + uE(u) + E(u 2 ) + E(v) + (E 2 (v))}- 
dt 3 

v[E~ l (v) + E-\u 2 ) + uE-\u) + u 2 + v}. (124) 

This is a Toda hierarchy on discrete time scales. The recursion relation 
between the n + 1-th and n-th elements of the hierarchy is given by 



v n+1 = uv n + vu n + vE-\u n )+v(E-\u)-u)(l- E)- 1 ^, (125) 
u n+l = E(v n )+uu n + v(l-E)- 1 — -E(v)(l-E)- 1 E—. (126) 

V V 

From this recursion relation the recursion operator of the hierarchy follows 

Example 7. Four-field system on time scale. We give two examples 
which are studied in [3] for the case T = Z. 

1. Let k = and a = —2 and 

L = £ 2 + w£ + v + u£~ 1 +p£- 2 . (127) 
Then we get the four-filed equations 

£=1 = vp-pE- 2 (v), (128) 

da 

— = vu + wE(p)-pE- 2 (w)-uE- l (v), (129) 
at i 

— = wE(u) + E 2 (p)-uE~ 1 {w)-p, (130) 
dti 

it = E >^- U - ( 131 ) 

2. Let k — 1 and a = — 2 and 

L = g£ 2 + w£ + t) + M£- 1 + £- 2 . (132) 
Then we get another four-filed equations 

nil 

i=\ = w-£~ 2 (w), (133) 
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wE(u) + q-E- 2 {q) -uE^iw), (134) 
wE(v) + qE 2 (u) - uE~\q) - vw, (135) 
qE 2 (v) - vq. (136) 

So far we considered the hierarchies coming from Proposition 19 with 
integer powers of the Lax operators. Now we consider the rational powers of 
the Lax operator. 

Proposition 22. Let 

L = w£ N + u N - 1 £ N ~ 1 + --- + u , (137) 

where w(x) is a function of x which is not a dynamical variable dw/dt = 0, 
Ui , i — 0, 1, • • • , N — 1 are functions of t and iGl Then 

^ = [(L"/ JV )> ,L], n = l,2,... (138) 

dt n 

produces hierarchies of integrable systems. Here n is a positive integer not 
divisible by N. Furthermore the function u is also not dynamical, i.e., uq = 
uo(x), not depending on t. 

Corollary 23 . When N = 2 and w = ^E(-) then the A-KdV Lax operator 
(81) reduces to the above form with 



dv 
dt\ 
dw 
~dT x 
dq 
~dT 



u = -- + \ + u, (139) 

ui = --[£(-) + -] + -. (140) 
fi fi fi fi 

Hence in part (2) of Example 5 we have a single equation with v = —fiu + 
l + fiu. 

In the following example we study the N = 2 case in more detail. 



23 



Example 8.KdV on discrete time scales. Let 

L = wE(w)£ 2 + u£ + v. (141) 

Then 

L 1/2 = w£ + a + a 1 £~ 1 + a 2 £' 2 + --- (142) 
where first three ctj are given as 

w(E(a ) + a ) = u, (143) 
toEfaij + E^Wai = v-{a ) 2 , (144) 

wE(a 2 ) + E~ 2 (w)a 2 = - ^"'^ ■ (145) 

Then we calculate L 3//2 by 

L 3/2 = wE(w)E 2 (w) £ 3 + p 2 £ 2 + pi£ +p + negative powers of £ (146) 
where 

p 2 = E(w)[wE 2 (a )+u], (147) 

Pi = wE(w) E 2 (a 1 ) +uE(a ) +wv, (148) 

p = wE(w)E 2 (a 2 ) +uE(a 1 ) +va , (149) 

7/ 7/ 

= wE-\w)[E- 1 + E{w)E}- 1 (E{a 1 )^)+v{l + E)- 1 ^-. (150) 

Then (138) with N = 2 produces a hierarchy of evolution equations. It turns 
out that v becomes a constant in the whole hierarchy. We give the first two 
members of the hierarchy (for n = 1 and n — 3): 



u tl = u(l -E)(l + E)~ 1 u ) (151) 
u t3 = u(l-E)p , (152) 

where po is given above. The next members of the hierarchy can be found 
by taking n = 5 in (138) or by applying the recursion operator TZ to u ts . For 
T = Kg and w — 1 the above hierarchy and its Hamilton formulation were 
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given by Frenkel [4]. The recursion operator of this hierarchy with w — 1 
can be found by using (58) with R n = a n £ + (5 n . We find that 



(E 2 -l)a n = E 2 (u n ), (153) 
(E 2 -l)(3 n = uE{u n ) + E{u)a n -uE{a n ) (154) 

and the equation which determines the recursion operator is 

u„+i = vu n - u(E - 1) /3 n , n = 0, 1,2, • • • . (155) 

We find that 

K = v - u(E + I)" 1 [-u + E(u)E](E 2 - l)' 1 E. (156) 

When the Lax operator is of degree one and has an infinite power series 
in operator the corresponding system is called the KP hierarchy. 

Proposition 24. Let 

L = S + w + wi£~ 1 +u 2 £~ 2 + ••• . (157) 

Then 

^ = [(L")> ,L], n=l,2,.-., (158) 
produces the following hierarchy 



n = l ^ = (E-l) Ul , (159) 



dt x 
dui 
~dt^ 

du k 
dt\ 



= (E-lfa + u^uo-E-^uo)], (160) 
= {E-l)u k+1 + u k [u -E~ k {u )], A; = 0,1, •••.(161) 



n = 2 ^ = (E 2 - l)u 2 -u 1 E~ 1 (E + l)u 



dt 



2 
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dui 
dt 2 



+E( Ul )(E(uo)+u ), 

(E 2 - 1) u 3 + a^uz) - u 2 E- 2 ( ai ) 

+a ui - u 2 E~ 1 (a ), 



(162) 



(163) 



du k 



(E 2 - l)u k+2 + a 1 E(u k+1 ) - u k+1 E k \a 1 ) 
+a u k - u k E~ k (a ), k = 0, 1, • • • , 



dt 2 



(164) 



where ao = (E + l)u\ + (wo) 2 an d cxi — (E + 1)m . The case T = Z of 
this hierarchy is discussed in [3] (see also the references therein) and the case 
T = K g is discussed in [4] and [5] . 

6. Trace Functional and Conservation Laws 

Let T be a regular time scale and A be the algebra of pseudo delta differential 
operators. Any operator F G A of order k has the form 

F = a k 5 k + a fc _i 5 k ~ l + • • • + ai 5 + a + a_i 5' 1 + a_ 2 <T 2 + • • • (165) 

where a/s are A-smooth functions of x G T (they are also functions of t G M). 
The coefficients ao and o_i we call respectively the /ree term (zero order term) 
and the residue of F associated with its " 5-expansion" (165) and write 



and therefore the same operator F can be expanded in series with respect to 
the powers of £ of the form 



F = b k £ k + 6 fc _x E h - X + --- + b 1 £ + b + b_ 1 £- 1 + 6_ 2 £~ 2 + ■ ■ ■ (168) 



Frees F = a (x) an d Res«5 F = a_i(x). 
In case of regular-discrete time scales T we have 



(166) 




(167) 
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We write 

Free^ F = b (x) and Res^ F = b-^x). (169) 
Substituting (167) and 

r 1 = (£- I)' 1 p = (S- 1 + £- 2 + -.-) f i = E-\fi) S- 1 + E- 2 (fi) £- 2 + ---, 

(170) 

into (165) and taking into account that 

E^ifi) = p{p{x)) = a(p(x)) - p(x) =x- p(x) = v(x), (171) 
we find that 

Res£ F — v Res<5 F. (172) 
Definition 25. The trace of an operator F G A is defined by 

Tr{F)= [ Res 5 {F(I + p5)- l }Vx, (173) 
where the nabla integral is defined according to Section 2. 

Proposition 26. Let F be given as in (165). In case of regular-discrete 
time scales we have 

Res 5 {F(I + pS)' 1 } = Free £ F (174) 

v{x) 

for x G T \ {x*,x*}, where x* = minT and x* = maxT. Therefore in this 
case 

( Free£F) Xl = J2b (x). (175) 

Proof. Since I + p5 = £ we have, by using (172) and (168), 

uRes s {F(I + P5)- 1 } = Res £ {F (I + P5)- 1 } = Res £ (F£~ v ) 
= Res^ (b k £ k ~ l + ••• + &! + b Q £~ l + •••) = b = Free £ (F). (176) 



Tr(F) = f 

JT 
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Proposition 27. For all F 1 ,F 2 E A 



Tr([F u F 2 \) = Tr(FiF 2 - F 2 F 1 ) = 0, (177) 
in other words the pairing (Fi,F 2 ) = Tr(FiF 2 ) is symmetric. 

We prove (177) only for particular cases of time scales T. They indicate a 
way to the proof in the general case of regular time scales, 
(i). If T = R, then 5 = d = £■ and pi(x) = 0, 

F = a k d k + --- + ai d + a + a^d' 1 + ■■■ (178) 

and 

Tr(F) = / Res{F)dx = / a^^dx. (179) 



It is well known that (for example, see [3]) for such functional Tr(F) the 
statement (177) holds. 

(ii). Let T be a regular-discrete time scale. Then by Proposition 26 we have 

Tt([F 1 ,F 2 ])= [ Res s {[F 1 ,F 2 }(I + f i5r 1 }Vx= [ (Free^, F 2 })^- = 0. 

Jt Jt v\x) 

(180) 

It is enough to check (180) for monomials F\ = AS k and F 2 = BS e . By the 
use of the property (113) of £ we have 

F X F 2 = A(E k B) S k+i and F 2 F X = B{E e A). (181) 

Therefore Frees[Fi,F 2 ] is either zero or 



Free f [Fi,F 2 ] = A{E k B) - B(E~ k A) = (E k - I)(E~ k A)B 

= (I- E- 1 ) {E k + E k ~ x + --- + E + I)(E~ k A)B = v(x)[$(A, B)f , 

where $(A, B) = (E k + E k ~ x + ... + E + I)(E~ k A)B. Hence 

/ Free £ [F 1 ,F 2 ]^- = [ [Q(A, B)]*Vx = B)£ = (182) 
Jt v \ x ) Jt 

so that (177) is proved for regular-discrete time scales. 
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(iii). Let T be a mixed time scale, say, of the form T = (—00, 0) U K 3 , where 
(—oo,0) denotes the real line interval. Then for any Fi,F 2 G A we have, 
taking into account Proposition 26, that 

Tr([F 1? F 2 ])= / Res 5 {[F 1 ,F 2 ](/ + /^)- 1 }V:r 
Jt 

= [ (Res d ([F 1 ,F 2 ])dx+ [ (Free £ [F 1 ,F 2 })^-. (183) 

J-00 JK q 

Take for instance Fi = A5 and F 2 = B5~ l . Then 



Res 9 [Fi, F 2 ] = Res d [Ad, BO' 1 ] = AB' + A'B = (AB)', 

Freest, F 2 \ = Free^ [- (£ - I), B(S - I)' 1 12} 

— AE(B) — BE~ 1 (A) (184) 



Therefore 



/0 rO 
(RjBB a [F 1 ,F 2 ])dx= / (AB)'dx = A(0-)B(0-), (185) 
-00 J —CO 

= [ A ( X ) B (Q X ) - A(q- 1 x)B(x)} = -A(0 + )B(0 + ). (186) 

X£K q 

Hence 

Tr([Fi, F 2 ]) = A(p-)B(p~) - A(0 + )B(0 + ) = 0, (187) 

where A and B are A-smooth functions on T and hence they are continuous 
at x — 0. 

Proposition 28 Equation (79) implies that 

^L k = [A n ,L% A n = (L%)> , (188) 
for all k = jj, where t is any positive integer. 
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Propositions 27 and 28 imply the next Proposition 
Proposition 29. For all i = 0, 1, • • • the junctionals 

Hi = Tr(L^), (189) 
are common constants of motion for the hierarchy (79) and (115). 

Note that in proof of Proposition 29 it is, in particular, used the fact that 
the flows (vector fields) defined by the different members of the hierarchy all 
commute with each other (see [2], [17]). 

8. Conclusion 

We have developed the Gel'fand-Dikii approach to time scales. So far the 
integrable systems were studied on R, Z or on K q . Here we gave a unified and 
extended approach. In particular cases when T = R, Z, K q we find several 
examples of the integrable systems. We developed the algebra of A-pseudo 
differential and £-shift operators. We established the GD formalism on these 
algebras and introduced several Lax representations on these algebras. All 
these Lax representations are straightforward generalizations of the Lax rep- 
resentations on pseudo differential algebras of integrable systems on R and 
the Lax representations of the algebra of shift operators on Z. The Burgers 
and KdV hierarchies on time scales that we found are the special cases of 
these Lax representations. We also generalized the Frenkel KdV hierarchy 
introduced on K g to arbitrary discrete time scales. We constructed the re- 
cursion operators of each example considered in this paper and gave a way 
of constructing the constants of motions by introducing an appropriate trace 
form on time scales. 

In this work we did not consider the r-Matrix construction and the Hamil- 
tonian formulation of integrable systems on time scales. The trace form on 
a general time scale needs a little care. Such a work is in progress and will 
be communicated in a separate paper. 
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APPENDIX: Recursion operators of four field 

systems 

We give the recursion operator of the Four-field systems on time scale which 
are studied in example 7. 

1. For the case k = 0, a = —2, we obtain the recursion relation between the 
n + 1-th and n-th elements of the hierarchy as follows: 

w n+l = w(E + l)- l E(v n ) + E(v)(E 2 - l)- 1 E 2 (w n ) - v(E 2 - 1)"V> + 

w(E + - E)w(E 2 - l)- 1 ^) + E 2 (u n ) + (1 - E 2 ) Vn ,(190) 

v n+l = wE(u n ) + vv n - u(E 2 - l^E-^Wn) + E(u)(E 2 - l)- 1 E 2 (w n ) + 

(1 - E 2 )p(l - E 2 )' 1 ^ + E 2 ( Pn ) + (E~ 1 (w) - wE) Vn , (191) 
u n+1 = wE{p n ) + u{E + l)" 1 ^ + 1 + E~ l )v n - p{E 2 - ly'E^iwn) + 
u(E + l)- 1 ^ - l)^- 1 ^)^ 2 - ly'wn + E{p){E 2 - l)' 1 E 2 {w n ) + 

(E- 2 (w) - wE)p(l - E 2 )- 1 ^ + vu n + (E-\v) - v)r] n , (192) 
Pn+1 = uE-\u n ) + p(l + E- 2 )v n + p{E- 1 - E~ 2 )w{E 2 - l)~ l E{w n ) + 

(E- 2 {v) - v)p{\ - E 2 )- lp -f + v Pn + {E-\u) - uE- 1 )^, (193) 

where 

Vn = (E 2 (p) - E{p)E 2 )^[E 2 {u)E{p n ) + E 2 {p)u n + {uE 2 {p) - E 2 {u)E{p)) 

(l-E 2 )-'E 2 (f)]. 

2. For the case k — 1, a — —2, the recursion relation between the n + 1-th 
and n-th elements of the hierarchy is given by 



u n+l = E~ 2 {w n ) + u(l + E)~\E - l)u(l - E^Eiun) + vu n + u(l + E)~ 1 v n + 

(E~\v) -v)(l- E 2 )-'u n + (1 - E~ 2 )( n , (194) 
w n+1 = uE~\q n ) + w(l + E)~\E 2 + E + l)v n + (E'^q) - qE 4 )(l - E 2 )~ l u n + 
w{\ + E)-\l - E)E{u){\ - E 2 )-^ 2 ^) + vw n + (E(v) - v)( n + 

(E 2 (u) - uE-^qiE 2 - l)- 1 ^! ), (195) 
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q n+1 = vq n + wE(w n ) + q(l + E 2 )v n + (E 2 (v) - v)q(E 2 - l)" 1 *? 2 ^ ) + 

q(l + E)-\l - E 2 )E(u)(l - E 2 )~ l E 2 {u n ) + (E(w) - wE)( n , (196) 
v n+ i = E~ 2 (q n ) + uE~\w n ) + (E~\w) - wE 3 )(l - E 2 y x u n + vv n + 

(E(u) - uE- l )( n + (1 - E 2 )q(E 2 - l^E 2 ^)., (197) 

where 

Cn = (^ 2 -^))-M^ 2 (^n) + (^ 2 (^)^-^(?))(^ 2 -l)- 1 ^(y)]. (198) 

From the recursion relations obtained in both cases, we can construct the 
recursion operators of the hierarchies. 



References 

[1] A.S. Fokas and I.M. Gel'fand, Bi-Hamiltonian structures and integra- 
bility, in Recent developments in soliton theory, edited by A.S. 
Fokas and V.E. Zakharov, Springer- Verlag, New York , 1992. 

[2] L.A. Dickey, Soliton Equations and Hamiltonian Systems, World 
Scientific, Singapore (1991). 

[3] M. Blaszak, Multi-Hamiltonian Theory of Dynamical Systems, 

(Texts and Monographs in Physics) Springer, Berlin (1998) 

[4] E. Frenkel, Internat. Math. Res. Notices, No.2, 55-76 (1996) 
(arXiv:q-alg/9511003). 

[5] B. Khesin, V. Lyubashenko and C. Roger, Journal of Functional Anal- 
ysis, 143, 55-97 (1997). 

[6] M. Adler, E. Horozov and P. van Moerbeke, Phys. Lett. A 242, 139- 
151 (1998). 

[7] B. Aulbach and S. Hilger, Linear Dynamic Process with Inhomogeneous 
Time Scale, in: Nonlinear Dynamics and Quantum Dynamical 
Systems (Gaussing, 1990), Math. Res., 59, Akademie Verlag, Berlin, 
1990, pp.9-20. 



33 



[8] S. Hilger, Results Math., 18, 18-56 (1990). 

[9] V. Kac and P. Cheung, Quantum Calculus, Springer Verlag, New 
York Berlin Heidelberg (2002). 

[10] M. Bohner and A. Peterson, Dynamic Equations on Time Scales: 
An introduction with Applications, Birkhauser, Boston (2001). 

[11] M. Bohner and A. Peterson, Editors, Advances in Dynamic Equa- 
tions on Time Scales, Birkhauser, Boston (2003). 

[12] F. M. Atici and G. Sh. Guseinov, J. Comput. Appl.Math., 141, 75-99 
(2002). 

[13] G. Sh. Guseinov, J. Math. Anal. Appl., 285, 107-127 (2003). 

[14] M. Giirses, A. Karasu and V.V. Sokolov, J. Math. Phys, 40, 6473-6490 
(1999). 

[15] M. Blaszak, Rep. Math. Phys, 48, no. 1-2., 27-38 (2001) 

[16] M. Giirses and K. Zheltukhin, J. Math. Phys. 42, 1309-1325 (2001). 

[17] G. Wilson, Math. Proc. Camb. Phil. Soc., 86, 131-143 (1979). 



34 



